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Abstract 



We consider the time dependent two point function, < <p q {t)<p_ q {0) > 

in non-linear stochastic field theories, for which the KPZ equation serves as a 
prototype, in particular we consider small q's and long times such that 0) q t » 1 (co q 

being the corresponding decay rate). We find that, since the generic case has 
(O q c< q^ 1 for small q where ji > 1 , the decay of the two point function is given by a 

stretched exponential in (O q t multiplied by a power of t, 

< <j)_ q (0)0,(0 ><* t Pd exp[-y(Q) q t)^] , where (5 d = (d - l)/2/j, , d is the 
dimensionality of space and y a dimensionless constant. 
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It is well known that many relaxation processes in condensed matter systems are 
characterized by stretched exponential decay. The list of examples contains 
relaxation in glasses 1 " 3 and in polymer systems 4 ' 5 , dielectric and viscoelastic 
relaxation ' etc. The theoretical derivation of stretched exponential decay involves 
usually extensive numerical simulations. This is true not only in cases where a 
realistic model is chosen to describe the physical system 8 but also in cases where 
more simplified models are used 10 " 13 . 

The present letter has a number of goals. The first is to show that phenomena of 
stretched exponential decay are much more abundant than can be expected from the 
canonical list presented above. Many interesting macroscopic systems such as 
magnets at the critical points, a turbulent fluid or a surface evolving under ballistic 
deposition etc., are non linear and share the property, that a disturbance of wave 
vector q decays with a characteristic decay rate 0) q , that scales at small q s as some 

power of q,0) q q* 1 . (Clearly, such systems attracted a lot of interest over a long 

period of time but that interest has not focused on the question of long time decay, 
that we address here). In most cases fi > 1 . An important exception is the case of 

turbulence where ji is about 2/3. We will show that for such systems that are 
generically described in terms of appropriate stochastic non-linear field theories and 
have ji > 1 , the decay is characterized by a stretched exponential function. The 
special case of turbulence will be discussed in the future. Our second goal is to 
construct a simple analytic framework within which the stretched exponential decay is 
derived. Our third goal is to show that the derivation is generic and that when the 
long time behaviour is considered most of the detail of the specific theory is not 
relevant. The only relevant parameter will be shown to be the exponent /i 

characterizing the small q dependence of the q dependent decay rate 0) q q* 1 . 

We consider a scalar field h(r,t) with Fourier transform h q (t) obeying a non-linear 
stochastic equation of the form 

dh 

"af +V A +L M ^>A« +LN qlmn h l h m K +■■=%, (i) 

where M,N have £ + m = q,£ + m + n = q etc. and f] q (t) is a noise term obeying 
<77 >=0 and <r](t)r)(t)>=2D°{q)8{t-?) (2) 



The KPZ equation, for example, has N=0 and M = 8 f+m _ q 7 ■ m—j= , where Q. is the 



volume of the system (to be taken eventually to infinity) and g is the non-linear 

coupling constant. For the dynamical <f> theory M=0 and N = S t+m+n _ q — (with the 

appropriate upper cut-off restriction) etc. The generalization of the above to vector 
fields needed in turbulence is straightforward. In order to be specific we will consider 
in the following the case where only M is different from zero but the derivation will 
readily suggest the treatment appropriate in other cases. 
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In our previous work 1415 we considered such systems and obtained the steady state 
"structure factor" ,<j> =< h q h_ q > s oc q~ r for small q and the corresponding decay rate, 

(O q , given by 



\<h q {t)h_ q {Q)>dt 

ft)" 1 = ^ — q -» for small q , (3) 



The method we employed is a self consistent expansion around a linear model that 
produces the "exact" tp q and co q . 

The expansion has thus the form 14 " 16 



(/> q =<p q +c q {(j) p ,(a p } , (4) 
explicitly, in second order of the expansion, 

, , 1 r r. / \ i V M tm (f) t (]) m M t M . 

<P q =<P q + —{D (q)-V q <p q + y + 2> q - , (5) 

ft), tf, 0) q + 0) f + 0) m ^ 0) q + 0), + 0) m 

and a similar relation for co 

0) q =0) q +d q {$ p ,a) p } . (6) 

The structure factor and decay rate are obtained by solving the coupled non-linear 
integral equations c {<j> ,co p } = and d q {(j) p ,a) p } = . In contrast to other 

expansions, the full correction, in a given order of the expansion, of the relevant 
physical quantities, is really small. In fact, it is chosen to be zero. In the following 
we obtain the "dynamical structure factor" <t> q (t) =< h q (0)h_ q (t) > s , using the same 

idea of a self consistent expansion. The average < — > s denotes that h q is measured 

in steady state at time t = and then h is measured at some later time t . The 

"dynamical structure factor" <J> ? (0 normalized by ? =<I> 9 (O), the static structure 

factor, is thus a measure of the persistence in steady state of disturbances with wave 
vector q . 

We will proceed first to obtain an equation for <J> ? (0 and then to show that in the 
small q and long time regime (co q t » 1) it must show a stretched exponential decay. 
Note that as it will become obvious from the derivation, the long time decay depends 
on the specifics of the problem only through the dependence of (i (co q ) on those 

specifics (as long as \x > 1 ). We derive the equation for <J> 9 (0 in detail but indeed the 
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most important part is the final simplified equation (11) that holds in the longtime 
regime. 



Our starting point is the field equation for h a , the Fourier transform in time 
h q (t) . 



+v q h qm + C qtm h fa h mt = r] qa , (7) 

where C qlm = — > T being an assumed periodicity in time to be taken eventually 

to infinity, a +T = 0) and the noise correlations are <r\ q J]_ q _ a >= 2D°(q) . (Note, 
that in the d+1 dimensional space that includes time, the noise is quenched disorder !) 

In the Chap man- Enskog spirit (as done in ref. 16) the equation is written in the form 
[(ifi) + co q )h q0) - rf qa ]+ X Yc qfm h f(S h mz -i$ +^ 2 [( v 9 - ft) 9 )V] =0 » ( g ) 

_£<T,mT 

where X is going to be taken as 1 but is used at present as an indicator to show the 
construction of the perturbation expansion as an expansion in X . The noise is split 

into two terms 77^ and 77^ such that <T] ( q °^ 77 _ ( °] m >= D qm and the correct <t> q(0 is 
given by <E> = — — T . This choice implies that ignoring the X and X 2 terms in 

9 CO :, I CO 

eq. (8) , we still obtain from a linear equation the correct <t> q(0 . Eq. (8) enables to 

obtain h qm explicitly to second order in X . The expression for h qco is multiplied into 

its complex conjugate and only terms up to second order in X are retained and then X 
is set to be 1 and the expression is averaged over i] qa) . The result is an equation of the 

form = ® qC o + e{® P o ) ■ Equating e{® pas } to zero yields 



W+v 2 q p qco -^c qem c eqm V— 4 J + r *— 

-2£lC 9 , ffl l 2 *^ fflT =D°( <? ) 



q(0 

(9) 



A _ 8 

Defining C qfm = — l ' +mq ; letting Q. and T tend to infinity and Fourier 

V £2 



transforming in 0) , we obtain 
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d 2 2 

V 

dt 2 " 



(2w) J J o 



-* q (t-n+a> q * q (t-n 



+ ]dt'e°" r ^ q _ l (f) 



-* q (t-n+a> q * q (t-n 



} 



^\d A i\A lA _ t \ 2 ® t (t)<t> q _ t (t) = D°(q)S(t) 



(10) 



1 r 

(Note that we define <£ (?) = — a e m d(0 and that <J> At) is symmetric in time). 

2k j 

The first step in the simplification of eq. (10) is to bring out Q q (t)of the integration 
over f . If we assume that 4> (0 falls off as a function of t slower than an 

exponential, it is possible to show that the most dominant contribution to the integral 
comes from the region with t' very small compared to t. This enables a 
simplification that leads indeed to the result that ^> q (t) decays slower than an 

exponential. Therefore the original assumption is consistent and the simplified 
equation reads 



a q ^ q {t)--^-^\d d i\A Lq _ f \% e (t)4> q _ e (t) = 



(11) 



where 



8ft) 



(27t) 



^\d d iA Lq _ e A qJ _ q \e-^ q _ e (t')dt' 



(12) 



and where we have neglected —r<I> (0 compared to <J>„(0, that is also consistent 

dt 2 q q 

with & q {t) decaying slower than an exponential. 

It is important to note here that there are many ways of obtaining equations for 
Q q {t) and indeed using different approximation schemes we obtained equations that 

differ from eq. (10). The important point however, is that their simplified form is 

similar to eq. (1 1) and leads to the same form of decay. 

The long time form of the solution is expected to be given by 



< f ? W = < ) / 1 (^)+<7 2 (®/) + 



(2) 



(13) 
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where lim = 



X— >oo 



/„(*) 



(Therefore, there is no reason to assume apriori thaty/ 1 ' =0 . In fact, if the series 
described by eq. (14) converges for all times (j) q =^V / ? ) if the f i 's are chosen to 

i 

have f t (0) = 1). 

The leading order form, f x , will be obtained now by balancing in eq. (11) the terms 
that are most dominant at long times. It is most important that apart from assuming 
that the appropriate integrals converge, the following discussion does not depend at 
all on the specific form of the M qfm (or equivalently A, A _ t ). 

Consider next a simple algebraic fact. The function A e {£ , q) = £ e + I q - £ f -q 
is negative in a small region of £ (for fixed q) foxO > 1 . (For <1 A e (£,q) is 
always positive). This implies that due to the fact that fi > 1 , /, cannot decay as an 
exponential or faster. If that were the case, the last term on the left hand side of eq. 

(11) would have a contribution from that small region of £ discussed above, that 
cannot be balanced against other terms in the equation, because it decays slower. 

To obtain the leading decay consider first the one dimensional case. The only 
possible way, by which the leading behaviour of the two terms left in equation (11) 

y 

can be balanced is by having f x (x) = exp[-]a /; " ] . The reason is that, say, for 

y y y 

positive q there is a finite segment < £ < q on which (O^ = + CO}/ 1 . (It is 
obvious that (O q °< g M only for small q but on the segment < £ < q small q implies 
small £ as well as small q - £ , therefore, (O, £/ and (O q _, (q- £Y ). 



The balancing condition leaves us with an equation relating the coefficients y/ 



{2k) j o 

The full solution involves the calculation of the coefficients y/ q l) from the equation 

above. The important point, however, is the stretched exponential form obtained for 
/, . The £ > q and £ < regions of integration produce a correction term of the form 

y 

exp[-ygr M ] times a negative power of t. The correction terms have to be balanced 
against correction terms to the time dependent structure factor (eq. (13)). 

It is immediately clear that a stretched exponential cannot solve eq. (11) for 
dimensions larger than 1 . If we take f (x) = exp[-jx s ] with 8 > ^/ the <J> 9 (0 
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term will dominate. If 8 < y the integral term will dominate, so that we are forced 

again at having 8 = ^j/ . The problem is that the relation [(0^ = [a> t ]^ + [&l ? 

holds only for vectors 1 of size bounded from above by q and parallel to vector q 

i.e., on a set of measure zero. Therefore, even for 8 = y,, , &„{t) dominates over 

/ Jl q 

the integral. The solution must, therefore, be of the form 

f l (x) = x p exp[-yx^] with /3>0. (15) 

To determine f5 , we have to consider the effective region of 1 integration in eq. 
(11). Let, 1 = l n +l ± , where l n is the part parallel to q and 2 ± the part 
perpendicular to it. The region of integration is of order q,t independent, in the 
direction of l n . The integration in the 1 ± directions is cut-off effectively at 

q 



■ oc t />l . This implies that the exponent [3 is given by 



■ (16) 

The derivation presented above focuses on the case where the non linear part of the 
stochastic field equation has in it only a product of two fields (like in the KPZ 
equation). When more than two fields are involved (like in the case of critical 
dynamics in the 4 theory) the situation is the same. The first term on the left hand 
side of eq.(ll) is linear <3> 9 (Oand the second term contains an integral on 
£ 1 ,£ 2 ,...,£ k aad the product <& n (t)<& e2 (t)...<& ek (t)<i> _ Vt (t) . The algebraic fact 

mentioned above is true for any number of variables. Namely, for a set of vector l i 

N 

obeying £4 =q, 

1=1 



£14 f -q e >0 for 9<\ (17) 

and for Q > 1 there is a small region of the £' s , for fixed q , where this expression 

becomes negative. Therefore, the derivation of the stretched exponential behaviour 
follows in the same way as above. The power (3 is also unaffected. It is easy to see 

using the arguments given above for the case of two fields in the non-linear part, that 
for n fields the equation determining is 

j8=/ij8 -4^(1-1). (18) 
2/j. 
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The long time form of <I> (0 is given therefore by 



® q (t) = c<l> q [nt / »] 2 exp[-r^], (19) 

where cis a dimensionless constant not necessarily 1!. The above form seems to be 
robust against various changes in the self consistent expansion and even against 
inclusion of higher orders. This will be discussed in future publications . 
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